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We show that the nonperturbative quantum transport 
equations, the so called 'Kadanoff-Baym equations', within 
the non-equilibrium real time Green's function description 
can be be understood as the ensemble average over stochastic 
equations of Langevin type. For this we couple a free scalar 
boson quantum field to an environmental heat bath with some 
given temperature T. The inherent presence of noise and dis- 
sipation related by the fluctuation-dissipation theorem guar- 
antees that the modes or particles become thermally pop- 
ulated on average in the long-time limit. This interpreta- 
tion leads to a more intuitive physical picture of the process 
of thermalization and of the interpretation of the Kadanoff- 
Baym equations. One also immediately understands that the 
emerging wave equations for long wavelength modes with mo- 
menta \k\ -C T behave nearly as classical fields. We also 
demonstrate how the problem of so called pinch singularities 
is resolved by a clear physical necessity of damping within the 
one-particle propagator. The occurrence of such ill-defined 
terms arising solely in a strictly perturbative expansion in 
out of equilibrium quantum field theory has a natural inter- 
pretation in analogy to Fermi's golden rule. 



I. MOTIVATION 

Non-equilibrium many body theory had been tradi- 
tionally a major topic of research for describing various 
transport processes in condensed matter physics and nu- 
clear physics. Over the last years a lot of interest for 
non-equilibrium quantum field theory has now emerged 
in particle physics. A very powerful diagrammatic tool is 
given by the 'Schwinger-Keldysh' [0 or 'closed time path' 
(CTP) technique by means of non-equilibrium Green's 
functions for describing a quantum system also beyond 
thermal equilibrium 0- The resulting causal and non- 
perturbative equations of motion (by various approxima- 
tions), the so called Kadanoff-Baym equations Q, have 
to be considered as an ensemble average over the initial 
density matrix characterizing the preparation of the ini- 
tial state of the system. Typically, if the system is close 
to thermal equilibrium, the (initial or resulting) density 
matrix allows for thermal fluctuations which should be 
inherent to the transport process under consideration. 
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If fluctuations are physically present in the course of 
the evolution, then, according to the famous fluctuation- 
dissipation theorem pj, also dissipation must be present. 
Or, in turn, if the (quantum) system behaves dissipa- 
tively, as a consequence, there must exist fluctuations. 
Again, the Kadanoff-Baym equations have to be under- 
stood as an ensemble average over all the possible fluctu- 
ations. This inherent stochastic aspect of the Kadanoff- 
Baym equations is what we want to point out and thus 
provide, as we believe, some new physical insight into its 
merely complex structure [§. Some of our conclusions 
are already known and can be found to some extent in 
various different studies [§Jf|. A new aspect, however, is 
the intimate connection of the Kadanoff-Baym equations 
to Langevin like processes. 

As a reminder of a Langevin process let us first briefly 
review the description of classical Brownian motion. 
Consider a classical system described by the variables 
Q and P interacting with a (heat) bath of oscillators of 
various frequencies pj according to the hamiltonian 
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From the resulting equations of motion the bath degrees 
of freedom can formally be integrated as 

q u (t) = q u (0) cos U) v t + - sin u> v t 

'rt toMt-m ^Qw (L1) 

o m v oj v 

so that the systems variables obey the effective Langevin- 
like equation of motion 

q + + 2 y o df r(t - t')Q(t') = m . (i.2) 

The bare hamiltonian H$ is slightly modified as 
H^\Q,P)=H s (Q,P)+(j2 



(1.3) 



On the other hand a Stokes-type frictional term (in gen- 
eral nonlocal in time) with a dissipation kernel 
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(1.4) 



enters into the left hand side of (1.2) supplemented by a 
source term £(i) on the right hand side being given as 
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£(t) = VrJ [g„(0) ^Q(O)] cosw** 

+ £^°l sill(J ,A . (1.5) 

In £ especially the initial conditions of the bath oscilla- 
tors enter. In general, they are not known exactly, but 
will follow some statistic distribution. Hence, has to 
be interpreted as a stochastic or 'noisy' source driving 
the fluctuations of the systems degrees of freedom. In- 
deed, if the initial conditions are randomly choosen by 
a thermal distribution P[q(0),p(Q)] ~ e - HB ^ BT §, 
is completely specified by a Gaussian distribution ('cen- 
tral limit theorem') with zero mean and the correlation 
kernel I 



«*(*)» =0, 



= 2k B TT(t - t') 



cosw„(i — t') 
(1.6) 



((...)) denotes the average over all possible realizations 
of the stochastic variable £(t). The dissipation kernel T 
and the random force £ (t) have the same microscopic ori- 
gin. In fact, as just derived, there exists a simple relation 
between them, called the fluctuation-dissipation theorem 
H . The evolution of the systems degrees of freedom are 
thus supplemented by an interplay between a dissipa- 
tive term and a fluctuating source. On the average, they 
should then equilibrate thermally. As a special appli- 
cation consider a heavy 'Brownian' particle with mass 
M placed in a thermal environment obeying an effective 
Langevin equation as just stated, i.e. 



Mx + 2 dt' T(t - 1') x(t') = £(t) . 



(1.7) 



(Here we assume that the friction kernel T(t) = T(—t) 
is 'well behaved' meaning that its Fourier transform 
T(lu) > and has only some finite extension in time.) 
To convince oneself that such a description in the long 
time limit is in accordance with the equipartition condi- 
tion {(p 2 ))/(2M) = T/2, one has to solve (Q by means 
of the retarded and advanced Green's function: 

GrctM - r— : p^y 



= g:». (i.8) 

One now finds in the long-time limit on the average the 
desired property 

«p 2 (^oo)» 



dhdt 2 G rct (t - ti)((Z(ti)S(t 2 )))G„(t2 - t) 
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G Tet W{u)G ay {w)=TM. 



(1.9) 



Hence, irrespective of the detailed form of the fric- 
tion kernel T(t), the equipartition condition is automat- 
ically fulfilled as long as the noi se k ernel lit) fulfills the 
fluctuation-dissipation relation (1.6). 

For a further physical motivation let us now return to 
quantum field theory and already point out some sim- 
ilarities. One of the major present topics in quantum 
field theory at finite temperature or near thermal equi- 
librium concerns the evolution and behavior of the long 
wavelength modes. These modes often lie entirely in the 
non-perturbative regime. Therefore solutions of the clas- 
sical field equations in Minkowski space have been widely 
used in recent years to describe long-distance proper- 
ties of quantum fields that require a non-perturbative 
analysis. A justification of the classical treatment of the 
long-distance dynamics of bosonic quantum fields at high 
temperature is based on the observation that the average 
thermal amplitude of low-momentum modes is large. For 
a weakly coupled quantum field the occupation number 
of a mode with wave vector p and frequency cop is given 
by the Bose distribution 



(1.10) 



For temperatures T much higher than the (dynamical) 
mass scale m* of the quantum field, the occupation num- 
ber for the low-momentum modes \p\ <C T becomes large 
and approaches the classical equipartition limit 
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> 1. 



(1.11) 



The classical field equations (which can be understood as 
a coherent state approximation to the full quantum field 
theory) should provide a good approximation for the dy- 
namics of such highly occupied modes. At a closer look, 
however, the cogency of this heuristic argument suffers 
considerably. The thermodynamics of a classical field is 
only defined if an ultraviolet cut-off k c is imposed on the 
momentum p such as a finite lattice spacing a. Many, 
if not most, thermodynamical properties of the classical 
field depend strongly on the value of the cut-off param- 
eter k c and diverge in the continuum limit (k c — > oo). 
In a correct semi-classical treatment of the soft modes 
the hard modes thus cannot be neglected, but it should 
incorporate their influence in a consistent way. In a re- 
cent paper of one of us || it was shown how to construct 
an effective semi-classical action for describing not only 
the classical behavior of the long wavelength modes be- 
low some given cutoff k c , but taking into account also 
perturbatively the interaction among the soft and hard 
modes. By integrating out the 'influence' of the hard 
modes on the two-loop level for standard (/) 4 -theory the 
emerging semi-classical equations of motion for the soft 
fields can be derived from an effective action and become 
stochastic equations of motion of Langevin type B : 
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(1.12) 



(Here (8> denotes a convolution in momentum space. The 
coefficients /i.; as well as the damping coefficients r)i de- 
pend on k and i, as well as on the momentum cut-off k c .) 
These resemble in its structure the analogous expression 
to eq. (1.7). The hard modes act as an environmental 
heat bath. They also guarantee that the soft modes be- 
come, on average, thermally populated with t he s ame 
temperature as the heat bath. Equivalently to ( |l.9| ) we 
will sec that a similar relation exists for the average fluc- 
tuation of the amplitude squared. For the semi-classical 
regime where \p\ <C T it will yield 
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(1.13) 



where V is the volume of the system. Such kind of 
Langevin description for the non-perturbative evolution 
of (super-)soft modes (on a scale of \p\ ~ g 2 T) in non 
Abelian gauge theories has recently been put forward || . 

In analogy to the Langevin description stated above we 
want to formulate in the following the effect of the heat 
bath on the evolution of the system degrees of freedom 
by means of the CTPGF technique. It is the intention 
of our study to provide in a self-contained way new in- 
sight into the dissipative and stochastic character of the 
CTPGF approach in thermal and non-equilibrium field 
theory. For this we discuss a free scalar field theory in- 
teracting with a heat bath. In the next section we first 
review the CTPGF technique and evaluate the average 
characteristic properties of the CTP propagator result- 
ing from the (non-equilibrium) equations of motion, the 
Kadanoff-Baym equations. Within a rather simple re- 
arrangement of the interaction kernels stemming from 
the heat bath, it becomes transparent to identify the dy- 
namically generated mass shift, the dissipation and the 
fluctuation terms, all of them contributing to the dy- 
namical evolution of the system. We also show that the 
fluctuation-dissipation theorem emerges naturally and is 
thus stated in microscopic terms. We then introduce in 
the third section the concept of a stochastic generating 
functional with external noise. In this way the analogy 
of thermalization in quantum field theory to a Langevin- 
like process becomes apparent. We thus can point out the 
influence of the noise correlation function on the propa- 
gator which contains the occupation number. As moti- 
vated above, one also immediately realizes that the long 
wavelength modes with momenta \k\ <C T and energies 
\u>\ <C T behave as classical propagating modes for a 
weakly interacting theory. We also briefly mention an at- 
tempt how to to derive an effective Boltzmann-Langevin 
equation for the phase-space occupation number of the 
bosonic particles from first principles. Section IV is de- 
voted to the issue of so called pinch singularities emerging 
in a strictly perturbative evaluation in out of equilibrium 



quantum field theory. We will give first a very clear phys- 
ical picture for their occurence. We then show how the 
'problem' of pinch singularities is cured naturally within 
the (non-perturbative) processes of thermalization and 
damping and providing a bridge to standard kinetic the- 
ory. We will end our findings with a brief summary and 
conclusions in section V. 



II. KADANOFF-BAYM EQUATIONS 

In this section we study the thermalization of a simple 
quantum field theory. For this purpose we couple a sys- 
tem of free scalar fields to an environmental heat bath 
of temperature T. We start with the closed time path 
action for this scalar field <fi ||: 



(-□ -TO 2 )c 

6+ £++ <h+ - 



- </T (-□ - m 2 ), 

>+£+- 4>- 

>-£— 6-} . 



(2.1) 



The time integration starts at some fixed time to at which 
the system starts to evolve from some initial density ma- 
trix. The interaction among the system and the heat 
bath is stated by an interaction kernel involving a self 
energy (or 'mass') operator £ resulting effectively from 
integrating out the heat bath degrees of freedom. Clearly, 
this self energy operator is the only quantity which might 
drive the system towards equilibrium. We assume that 
£ is solely determined by the properties of the heat 
bath. Ideally such a scenario holds for the linear response 
regime. 



In (2.1) the self energy contribution from the heat bath 
is parametrized in the Keldysh notation by the four self 
energy parts 



£ ++ (xi, x 2 ) = 

Q{t 1 -t 2 )Y 1 > {x ll x 2 ) + Q{t 2 ^t 1 )H < {x ll x 2 ) 
£ (xi,x 2 ) = 

e(t 2 -t 1 )Y, > (x 1 ,x 2 ) + e{t 1 -t 2 )T J < (x 1 ,x 2 ) 

Y> + ~(xi,x 2 ) = -£<(:£!, £2) 
Y>~ + (xi,x 2 ) = -£> (xi,x 2 ) 



(2.2) 

(2.3) 
(2.4) 
(2.5) 



where ti and t 2 are the time components of the vectors 
x\ and x 2 , respectively. With these definitions at hand 
one sees that in principle only two self energy parts are 
independent, e.g. £ < and £ > . One can also easily show 
that the difference £ > (a;i,a;2) — Y> < {x\,x 2 ) is real while 
the sum £ > {x\, x 2 ) + T, < (xi,x 2 ) is purely imaginary. 

If, as assumed, the heat bath stays at thermal equilib- 
rium at some fixed temperature T and the self energy is 
solely determined by the property of this heat bath, then 
the important relation 



£>(fc) = e fco/T £<(fc). 



(2.6) 
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holds. This property (2.6) which is nothing but the 
Kubo-Martin-Schwinger boundary condition flgJlJill 
will become crucial below when we will show that the 
modes of our system in the long-time limit will indeed 
equilibrate at the temperature of the heat bath. 

In order to explicitly explore the causal structure, it is 
useful to introduce the more physical retarded and ad- 
vanced self energies 

X Iet { Xl ,x 2 ) :=0(ti-t2) [^ > (x 1 ,x 2 )~^ < (x 1 ,x 2 )] , 
Z a - V (x u x 2 ) :=0(t 2 -*i) [^{xux^-^ixux^] . 

Given the fields <fi + and <\>~ on the two branches of the 
CTP contour, four two-point functions can be defined: 

D ab ( Xl ,x 2 ) := -i{r(xi)^ b (x 2 )} , 0,6=+,-. (2.7) 

The expectation value is defined via path integrals as 



(-□ 



D< - E rct D< - E<£) a 



0, 



(2-15) 



(O) := i /p^^lOe^^^M^.f] (2.8) 

where p is the density matrix of the system for the initial 
time to. 



When evaluating the two-point functions by using (2.6 
one finds 



D++{x u x 2 ) = 

0(ti - t 2 ) D > (x 1 ,x 2 ) + Q(t 2 - h) D<( Xl ,x 2 ) , (2.9) 
D (2:1,2:2) = 

0(t 2 - h) D > {x 1 ,x 2 ) + 0(ti - t 2 ) D<( Xl ,x 2 ) , (2.10) 
D+-{ Xl ,x 2 ) = -i{<p-{x 2 )<P + ( Xl )) =: D<{x u x 2 )i2.11) 
D-+{ Xl ,x 2 ) = -i{<p-{ Xl )^+(x 2 )) =: D>( Xl ,x 2 )(2.12) 

In addition, one also now introduces retarded and ad- 
vanced quantities by 

D let (x 1 ,x 2 )=e(t 1 -t 2 ) [D>(x 1 ,x 2 )-D < (x 1 ,x 2 )] , 

£ av (2:1,2:2) =e(t 2 -ti) [r>< (2:1, 2:2) -r>> (11,12)] . 



From (2J) it is straightforward to obtain the equations 
of motion for the matrix elements. Since the four two- 
point functions are not independent of each other the four 
equations will contain redundant information. Therefore 
we will state the equations of motion for retarded and 
advanced propagator and for an off-diagonal element, say 



D< := D 



(2.13) 



The equation of motion for the retarded propagator then 
reads 



-□ 



)D* 



(2.14) 



A similar one is given for the advanced propagator. 
Please note already the analogy to the solving of the 
Langevin equation carried out in the introduction by in- 
troducing the retarded and advanced propagator. Addi- 
tional information now comes from the equation of mo- 
tion of the propagator D < , which reads 



Th is is j ust th e fam ous Kadanoff-Baym equation ||. 

( 2.14 ) and (|2.15| ) determine the complete and causal 
(non-equilibrium) evolution for the two-point functions. 
To get now more physical insight into the (effective) ac- 
tion given in ( |2 . 1[ ) and in the resulting Kadanoff-Baym 
equations we now introduce the following real valued 
quantities: 



s{x\,x 2 ) 



a{x\,x 2 ) 



I(xi,x 2 ) 



t 2 ) (E > (2;i,2;2) 



1 , 

s(x 2 ,xi) , 
- {TP > {x 1) x 2 ) 
-a(x 2 ,xi) , 

-— (E > (xi, x 2 ) + (xi, x 2 )) 
I(x 2 ,xi) . 



S < (2;i, x 2 )) 

(2.16) 



S < (a3i,x 2 )) 



(2.17) 



(2.18) 



Our notion for s and a serves as a reminder for the re- 
spective symmetry properties. It basically represents the 
standard decomposition of the real and imaginary part 
of the Fourier transform of the retarded self energy op- 
erator E rct s yields a (dynamical) mass shift for 
the <f> modes caused by the interaction with the modes of 
the heat bath, while a is responsible for the damping, i.e. 
dissipation of the <fi fields. The important thing to point 
out will be that / characterizes the fluct u ation s. 

We first note that on account o f (2.16 )-( 2.1S| ) together 
with (2. 2)- (2. 5) the CTP action (2.1) can be written as 



S 



1 



+ ii 



- (j) ) (s + a) 1 
+ -<P-)I{<P+ 



+<n 



(2.19) 



This expression is identical to the influence functional 
given by Feynman and Vernon JT^Jl3] |. To the exponen- 
tial factor in the path integral (|2.8| ) the (s + a) term con- 
tributes a phase while the I term causes an exponential 
damping and thus signals nonunitary evolution. Follow- 
ing the ideas of Feynman and Vernon one can therefore 
indeed identify I as a 'noise' correlator (see next section) . 
The equivalence of the CTP formalism and the influence 
functional approach has indeed been pointed out already 
in |0| on a strictly formal level. Our discussion here is 
inspired on more physical intuition. 

The two relevant equations of motion are stated as 

(-□-m 2 - s-a)D Ict =6, (2.20) 
(-□ - m 2 - s - a)D < + (a + H)D av = . (2.21) 

We see that the last equation (the 'Kadanoff-Baym' equa- 
tion) is the only one where I occurs. Hence the retarded 
(and advanced) propagators are determined by s and a 
while the number density is additionally influenced by 
the fluctuations /. 
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For the interpretation of s and a we consider the long- 
time behavior of these equations. In this case we can 
assume that the system becomes translational invariant 
in time and space and the boundary terms are no longer 
important. For £> ret one immediately finds 



It is now very illuminating to explicitly write down the 
relation between a(k) and I(k) 

^' Iw-s'w '^'^)"'"' (2 - 28) 



D iet (k) 



1 



fc 2 — to 2 — s(k) — a(k) 
From this the spectral function is given by 



D av (fe))* (2 22) ^ n ^ nc high temperature (classical) limit one gets 



A(k) 



-ID^(k) 



-D av {k)} 
i a(k) 



[k 2 - m 2 - s(k)] 2 + \a(k)\ 



(2.23) 



Note that a(k) is purely imaginary since a(x\, X2) is real 
and antisymmetric. For the same reason s is real. In- 
specting the spectral function it becomes obvious that 
s = ReS ret contributes an (energy dependent) mass shift 
while a = iImE ret causes the damping of propagating 
modes j3|(| . a is related to the commonly used damping 
rate T via 



f(fc) 



.a(fc) 
k 



(2.24) 



For D < one finds in the long-time limit the relation 

D<{k) = D Iet (k)t, < D sv (k) 

= D rct (k)[-a(k) - iI{k)]D m {k) 
-a{k) -il(k) 



A{k) ee -2in{k)A(k). (2.25) 



We note that the first equation stated in (2.25) is some- 
times denoted in the literature |||| as a generalized 
fluctuation-dissipation theorem. We will now outline its 
intimate connection to a standard Langevin process in 
the following. 

To proceed we need a relation between a and /, i.e. a 
relation between the damping and the noise term. In 
our case it is a simple cons equen ce of th e KMS condition 
( f2~6[) using the definitions (|2.17D , (|t|): 



n(k) 



£<(*) 



f 



E>(fc)-E<(fc) e k °/ T -l 



n B (k ), (2.26) 



which indeed shows that the phase space occupation 
number of the soft modes in the long-time limit becomes 
a Bose dist ribution with the temperature of the heat 
bath. ( [2.26] ) is independent of the explicit and detailed 
form of the self energy but is solely determined by the 
KMS condition. If one now assumes that the coupling 
between the bath and the sys tem becomes very weak, 
i.e. a, f , J — > 0, the expression (2.26|) 



n(k) 



'0' 



1 



e ko/T _ I 



(2.27) 



is still preserved as long as the KMS condition is fulfilled. 



I(k) = -^2ia(k) , 
k 



or, employing ( 2.24 ), 

I{k) = 2Tf(fc) 



(2.29) 



(2.30) 



Recalling our discussion of Brownian motion in the intro- 
duction this compares favorably well with (1.6)! Indeed 
as we shall see in the next section we can define a quan- 
tity which obeys a Langevin equation very similar to the 
one for the Brownian particle. The physical meaning of / 
as a ' noise' correlator will become obvious. The relation 
( [2.28 ) thus already represents the generalized fluctuation- 
dissipation relation |7J from a microscopic point of view 
by the various definitions of /, a and T through the parts 
E < and E > of the self energy. 

We close this section by remarking that in all appli- 
cations typically the major goal of the Kadanoff-Baym 
equations had been to derive a standard kinetic transport 
equation for the (semi-classical) phase-space distribution 
f(x,k,t) which should be valid for weak coupling and a 
nearly homogeneous system. If the system is in a general 
off-equilibrium state the two-point functions depend not 
only on the so called microscopic variable u :— x^—x^ but 
additionally on the macroscopic center-of-mass variable 
X := (xi + X2)/2. Performing a standard Wigner trans- 
formation for D K and making a gradient expansion up to 
linear order in the center-of-mass variable X one finds the 
well-known relativistic kinetic transport equation in the 
quasi-particle approximation for the semi-classical (on- 
shell) one-particle phase-space distribution HJ|,|| 



kftd^fix, k,t) = 

~ (it< (X, k) [f(x, k, t) + 1] - zE> (X, k) f(x, k, t) 



(2.31) 



k =u°. 

resembling in its form the Boltzmann equation. Here all 
energies have to be evaluated onshell. In this form it is 
obvious that iE < /2w2 can be interpreted as the produc- 
tion rate, while zE > /2o;2 stands for the absorbtion rate 

' k 

for modes with the respective energy. 



III. STOCHASTIC INTERPRETATION AND 
SOME CONSEQUENCES 

To see the connection between the formalism presented 
in the previous section [n] and stochastic equa tions we 
decompose the influence action S as given in ( [2.19 ) in 
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its real and imaginary part and rewrite the generating 
functional 

z\j + ,r] 

x e iReS[<l>+,<l>-]+ij+<t>++ij-<t>--i(4>+-$-)I(4>+-$-) 

= i Jme-^ 1 '^ J-D[4> + ,4>-]p[<f> + ,4>-] 



x e 



I f-D[Z]e-K I - 1 *Z'\j+ + U--Z] 



(3.1) 



with 



N := V£e 



The generating functional Z[j + ,j ] can thus be inter- 
preted as a new stochastic generating functional Z'[j + + 
— £] averaged over a random (noise) field £ which is 
Gaussian distributed with the width function 7, i.e. 



«0» 



iv7 



(3.3) 



From the last definition we find that the (ensemble) av- 
erage over the noise field vanishes 



= 0, 



while the noise correlator is given by 



(3.4) 



U0) = i- 



(3.5) 



The action entering the definition of Z' is no long er S, 
but only the real part of the influence action ( 2.1 9| ) . 

From this new stochastic functional Z' a Langevin 
equation for a classical <j> field can be derived ||. Not- 
ing that the fields (</> + ){ on the upper branch and (<f>~)^ 
on the lower branch are equal (and denoted as 4>£, m the 
following), its equation of motion derived from Z' takes 
the form 



(-□ - m 2 - s) -afe = 



(3.6) 



This, indeed, represents a standard Langevin equation. 
The spatial Fourier transform of the Langevin equation 
(|3.6[) then takes the form 



fa(k,t) 



2T(k,At = Q))<p 6 (k,t) 



+2 / dt'T(k,t-t')^{k,t') =^(k,t), (3.7) 



The analogy between this Langevin equation (3.7) and 
the one for a single classical oscillator is obvious. The 



important differenc e, how ever, is th e fact that the corre- 
sponding relations ( 2.28 ) and (1.6) between the respec- 
tive noise kernel I and friction kernel T only agree in the 
high temperature limit. 

One can now ask to what ext end the classical equations 
of motion (3J3) together with (3.5) are an approximation 
for the full q uantum problem giv en b y the equation of 
motion ( 2.21 ) for D < . Inverting (3.6) one finds for the 
correlation function in the long-time limit 

- i(( <0+) 5 (0-> { » = -i D* et (mD™ = -i D Iet ID av . (3 



Note that ( |3.8| ) is indeed the relation ( 1.1 3| ) advocated 
in the introduction to hold in the (semi-)classical regime. 
This one has to com pare with the full quantum correla- 
tion function D < of (2.25). One thus has that (—a — il) 
is approximated by —il. Of course this is justified, if 
|a| C J holds. Using the microscopic quantum version 



(3-2) ( pT28|) of the fluctuation-dissipation theorem this is equiv- 
alent to coth (|^) 3> I- Thus in the high temperature 
limit or - turning the argument around - for low fre- 
quency modes, i.e. for fco <C T, the classical solution 
yields a good approximation to the full quantum case. 
To be more precise: In simul atio ns one has to solve the 
classical Langevin equation ( |3.6| ) and calculate n-point 
functions by averaging over the random sources. This 
has been raised in Q. We would also like to refer to 
the works in |l5f| for some practical examples regarding 
the mathematical evaluation of expectation values in the 
transition from quantum to classical field theories using 
dimensional reduction techniques. 

One can also write down the equations of motion for 
the quantum two-point functions with external noise || 
by introducing the 'noisy' propagators 



Df{ Xl ,x 2 ) 



-z(0 o (a: 1 )/(x 2 ))£ 



(3.9) 



One recovers that £>^ et and D| v obey the same equations 
of motion as D ret and D av , respectively, and are thus the 
same. Only the equation of motion for and hence for 
the occupation number is modified compared to the one 
for D < . Av era ging this equation over the noise fields ac- 



equation (2.21 



cording to (3.2) one indeed rederives the Kadanoff-Baym 



This demonstrates that the Kadanoff- 
Baym equation can be interpreted as an ensemble aver- 
age over fluctuating fields which are subject to noise, the 
latter being correlated by the sum of self energies S < 
and E > , i.e. from a transport theoretical point of view 
the sum of production and annihilation rate. We want 
to note once more that the 'noisy' or fluctuating part 
denoted by I in heren t to the structure of the Kadanoff- 
Baym equation ( 2.15 ) guarantees that the modes or par- 
ticles become correctly (ther mall y) p opulated, as can be 
realized by inspecting ( 2.26 ) or ( 2.27 ). 

What is changed, if we replace our toy model of a free 
system coupled to an external heat bath by a self-coupled 
and thus nonlinear closed system? In an interacting field 
theory of a closed system the Kadanoff-Baym equations 
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formally have exactly the same structure as in our toy 
model. The important difference, however, is that the 
self energy operator is now described fully (within the ap- 
propriate approximative scheme) by the system variables, 
i.e. it is expressed as a convolution of various two-point 
functions. Hence, an underlying simple stochastic pro- 
cess, as in our case an external stochastic Gaussian pro- 
cess, cannot really be extracted. However, we emphasize 
again that the emerging structure of the Kadanoff-Baym 
equations is identical. The decomposition of the self en- 
ergy operator into its three physical parts (mass shift s, 
damping a, and fluctuation kernel 7) can immediately 
be taken over. Hence these three parts keep their clear 
physical meaning also for a nonlinear closed system. 

We close our discussion by noting that one can also 
pursue to derive a standard kinetic transport equation 
for the (semi-classical) phase-space distri butio n fix, k, t) 
including fluctuations H. In analogy to (2.31) one gets 



k^f^xXt) = ~ (it<(X,k) [f 6 (x,k,t) + 1] 



-it>(X,k)f s (x,k,t) 



kQ—LJ l 



+ u)jft(x > k,t). (3.10) 



This derived kinetic transport process ( 3.10 ) has the 
structure of the phenomenologically inspired Boltzmann- 
Langevin equation ||lq| : In order to describe fluctua- 
tions around the average, Bixon and Zwanzig Jl6[ | postu- 
lated a stochastic (classical) Boltzmann equation in anal- 
ogy to the Langevin equation for a Brownian particle. 
By adding a fluctuating collision term to the linearized 
Boltzmann equation (of the form of equation (2.31)) they 
obtained the following scheme 



dt 



fdt) = -nfdt)-.f cq )+F? hcn (t) 



(3.11) 



for a system near equilibrium. The correlation function 
of the fluctuating collision term was guessed on the basis 
of the fluctuation-dissipation theorem. Also it was shown 
that the stochastic Boltzmann equation provides a basis 
for describing hydrodynamic fluctuations. On the other 
hand, the equations of Bixon and Zwanzig were not de- 
rived but obtained on the basis of intuitive arguments. 
Instead, our approach carried out in Q has to be consid- 
ered as a clear derivation from first principles. Indeed it 
shows (nearly) a one to one correspondence to the phe- 
nomenologically introduced scheme. However, also some 
severe interpretational difficulties in the interpretation of 
the fluctuating phase-space density remain. We refer the 
interested reader to our discussion in Q . 



IV. NONPERTURBATIVE RESOLUTION OF 
PINCH SINGULARITIES 

We have seen in the previous sections that modes 
or quasi-particles become thermally populated by a 



non-perturbative interplay between noise and dissipative 
terms entering the Kadanoff-Baym equations. It is non- 
perturbative as the equations of motion explicitly resum 
the self energy contributions by a Schwinger-Dyson equa- 
tion defined on the real time path contour. In addition 
we have observed in eq. ( 2.27 ) that seemingly ill-defined 
expressions of the form '0'/'0' are well-behaved in the 
sense of a weak coupling limit. 

Rather similar ill-defined expressions result from so 
called pinch singularities within the context of strictly 
perturbative expansions in real time non-equilibrium 
field theory first raised in ref. jl7|. It is the purpose 
of this section to demonstrate explicitly how pinch sin- 
gularities are absent within a necessary non-perturbative 
context |Q. 



-e- 



FIG. 1. Lowest order self energy term in <^ 4 -theory which 
contributes to the pinch problem (sunset diagram) . The imag- 
inary part of the sunset diagram can be identified with a scat- 
tering amplitude. 

Employing the diagrammatic CTP rules potential 
'pinch singularities' can arise in strictly perturbative ex- 
pressions. Consider in the following a weakly interacting 
scalar </> 4 -theory. The initial state in the far past (assum- 
ing a homogeneous and stationary system) is prepared 
by specifying the momentum occupation number n(p) of 
the (initially non interacting) onshell particles. Note that 
this occupation number depends only on the three mo- 
mentum p. The occupation number n(p) enters the (free) 
propagator 

d o(p) = -27risgn(p ) 5{p 2 - to 2 ) 

x [e( Po )h(p) - e(- Po )(l + n{p))\ (4.1) 

The free retarded and advanced propagator are given by 



D 



ret / av 




(p) 



i 



p2 _ m 2 _|_ j £ S g n (p ) 



(4.2) 



A specified self energy insertion S p in a strictly pertur- 
bative expansion is given by a convolution of these free 
propagators. As a characteristic example the 'sunset' 
graph arising in scalar 4 -theory is illustrated in fig. [l]. 
We choose this particular graph as an example since the 
self energies Ep^^po = E p ) do not vanish onshell. 
The strictly perturbative corrections to the propagators 
within the context of the Schwinger-Keldysh formalism 
take the form 



£jrct/av _ £ji'ct/av _j_ jy et / av jjret / av £) re V av _|_ 



=:D 
D< = D< 



ret / av 



7-,rctv< n av 
u % U Q 



7->< v av n av 
U Q 



=:D< 



AD < 



(4.3) 
(4.4) 
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As £>Q Ct contains a pole at po = ±E p — it and Dq v a pole 
at po = ±-Bp + ie the product of both in the expression 
for AD< is ill-defined, if E<(p,p = E p = yjm? + p 2 ) 
does not vanish onshcll. Transforming such an expres- 
sion back into a time representation, the contour has to 
pass between this pair of two infinitely close poles. In 
fact, all three contributions to AD < are ill-defined. On 
the other hand the perturbative corrections A£) ret / av to 
the free retarded/advanced propagator are free of any 
pinch singularities as the emerging poles are all located 
at the same side of the contour. The three perturbative 
corrections can be further rearranged jl7]|| to extract the 
part carrying the pinch singularity 



A£>< inch (p) 



(4.5) 

D T et (p) [0(po) ((1 + n(p))E<(p) - n(p)E>(p)) 

+e(-po) ((1 + n(p))E>(p) - n(p)E<(p))] D^(p) . 

The expression in the square brackets is familiar from 
standard kinetic theory (compare with 2.31). Apart from 
a trivial factor one can interpret 



1 



r cff (p) := [(l+n(fl)»E>(p)-nGS)iE<. 



Po=E p 

(4.6) 



as the net effective rate for the change of the occupation 
number per time. For an equilibrium situation the oc- 
cupation number is given by the Bose distribution an d 
the self energy insertions fulfill the KMS condition (2.6). 
Hence, for the equilibrium case the whole bracket ex- 
actly vanishes and no pinch singularities emerge within 
a strictly perturbative expansion. It was observed and 
proven already by Landsman and van Weert that such 
ill-defined terms do cancel each other, to each order in 
perturbation theory, if the system stays at thermal equi- 
librium [n]. This is not the case for a general non- 
equilibrium configuration This seemingly 'severe' prob- 
lem arising for systems out of equilibrium was first raised 
by Altherr and Seibert . 

In a subsequent paper Altherr jl9| tried to 'cure' this 
problem by hand by introducing a finite width for the 'un- 
perturbed' free CTP propagator D so that the expres- 
sions are at least defined in a mathematical sense. Such 
a procedure, of course, already represents some mixing of 
non-perturbative effects within the 'free' propagator. As 
we have seen from our discussions in the previous sections 
a damping term (resulting in a non vanishing width) and 
the associated noise guarantee that the propagator be- 
comes thermally populated at long times. Hence, for any 
system which moves towards thermal equilibrium and 
thus behaves dissipatively, the full propagator must have 
some finite width (due to collisions or more generally due 
to damping) . Within his modified perturbative approach, 
Altherr |1!J also showed that seemingly higher order di- 
agrams do contribute to a lower order in the coupling 
constant, as some of the higher order diagrams involving 



pinch terms will receive factors of the form 1/T n , n > 1 
reducing substantially the overall power in the coupling 
constant. In his particular case Altherr investigated the 
dynamically generated effective mass (the 'tadpole' con- 
tribution) within standard 4 — theory. (For the hard 
modes the onshell damping T is of the order of o(g i T).) 
Therefore he concluded that power counting arguments 
might in fact be much less trivial for systems out of equi- 
librium. We will come back to his observation below. 



Comparing (4.5) and (4.6) with the Boltzmann equa- 
tion (2.31) one is already tempted to speculate that 
the occuring pinch singularity must be connected to the 
change in time in occupation number. Hence, we will first 
now elaborate on the physical reason for the occurrence 
of pinch singularities. Within standard scattering theory 
one would think about the probability for a particle of 
some initial momentum state to be scattered into another 
momentum state. Therefore we ask, how the occupation 
number n has changed after a long time. The occupation 
number for the out-states can be readily extracted from 
D < by means of the formula [n~8j 



n(p,t 



j (out) = (fl t (out) (out) ^ 



E„ 



1 d d 
2E D 'didt' 



d 3 ye l 



i d d 
2^dt ~ &F 



1 

V 



(iD<{y,t;x,t')) 



(4.7) 



We now assume that the interaction is switched on at 
a time t = —T/2 and switched off at t = T/2, i.e. we 
replace Ep ^ [x\, x%) by 

e(f - ti) e(f - 1 2 ) t.<i>{x u x 2 ) Q(h + f ) &(t 2 + f ) 

and assume that the duration time T is large but finite. 
This procedure regulates the pinch singularities to a finite 
value. Invoking a couple of manipulations p8| one ends 
with the rather 'obvious' result 



« r eff (p1 / d £ 
«r eff (p)-T. 



(4.8) 



We thus have demonstrated the bridge between the oc- 
currence of pinch singularities within the context of the 
CTP formalism and Fermi's golden rule in elementary 
quantum scattering theory. The effective rate F e ff is 
therefore analogous to the transition probability per unit 
time. Indeed one can easily understand in physical terms 
that one has to expect such a singularity in perturba- 
tion theory: Staying strictly within the first order contri- 
bution the particles remain populated with the initially 
prepared non-equilibrium occupation num ber (since this 
quantity enters the free propagator (hi)) and scatter 
for an infinitely long time. Therefore, the resulting 
shift Ari(p)(° ut ) should scale with T cfi (p) ■ T with T cS (p) 
held fixed. However, looking at the Boltzmann equa- 



tion (2.31), the occupation number and the collision rate 
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do not stay constant during the dynamical evolution 
of the system, but will be changed on a timescale of 
roughly 1/r. The quasi-particles are not really assymp- 
totic states. 

Subsequently, pinch singularities are formally cured by 
a resummation procedure. The propagators will become 
dressed and supplemented by a finite (collisional or more 
generally damping) width. This represents already a non- 
perturbative effect which only can be achieved by a re- 
summation of Dyson-Schwinger type. As a first attempt 
(proposed by Bai er et al. pO|), one might resum the full 



series of (13jT4) using the self energy E p . With the 



definitions T 
ReS 



P (p,Po) 



Re£ p ot 



ReEp V one 



finds @ 



D 



= (-2i) 



ret ^ D ^)(-^) D <( D ^)(-^ 
Po^p 



D a 



S p(P>Po)] and 



£> ret £<D av 



(p 2 



171 



(4.9) 



-ns(p>o) 



(The first term on the r.h.s. of the first equation rep- 
resents a boundary term |5||| of the (initial) propagator 
. If there would be no interaction, i.e. £ p = 0, then 
this boundary term just gives the free and undisturbed 
propagator. On the other hand, if S p is nonvanishing 
onshell, as assumed, this boundary term vanishes within 
the carried out resummation Only the second term 
then contri bute s and yields the result stated in the sec- 
ond line of (|4~|).) 



Hence the resummation of the series (4.4) of ill-defined 
terms results in a well-defined expression. The quantity 
n^(p, po) has to be interpreted as the 'occupation num- 
ber' demanded by the self energy parts |J. If the equilib- 
rium KMS conditions (2.6) apply for the self energy part, 
then ns(p,po) — > ns{po) becomes just the Bose distri- 
bution function. For a general non-equilibrium situation, 
however, thi s fa ctor de viat es from the Bose distribution. 
Comparing (11) with (^j|), the astonishing thing to note 
is that in fact the (initial) non-equilibrium distribution n 
has dropped out completely! 

Calculating £ p on a purely perturbative level the ini- 
tial occupation number h enters via the free propaga- 
tor (4.1). This however cannot be the whole truth in a 
dynamically evolving system. It is important to make 
sure that such a system is prepared at some finite initial 
time to- (If t would be taken as t — > — oo the system 
would already have reached equilibrium long time ago. 
Time reversal symmetry is explicitly broken, so that the 
propagators in principle have to depend on both time 
arguments explicitly before the system has reached a fi- 
nal equilibrium configuration. Therefore the simple use 
of Fourier transforms, as carried out in |l7Ul9lj20|l, and 



which in fact has led to the pinch singularities in (4.5)), 
is dubious |Ej.) The initial out of equilibrium distribu- 
tion n{to) cannot stay constant during the evolution of 
the system as it has to evolve towards the Bose distri- 
bution. As long as the system is not in equilibrium, the 



propagator thus cannot be stationary. In addition, the 
self energy parts E < and S > do also evolve with time. 
Hence they should depend on the evolving distribution 
function and not persistently on the initial one, ft, which 
enters S p in (4.£). Thus, the resummation of (4.9) docs 
not cover all relevant contributions. The self energy op- 
erators must be evaluated using the fully dressed and 
temporally evolving one-particle propagators. 




FIG. 2. Dyson-Schwinger equation with fully dressed prop- 
agators (skeleton expansion). 

The solution to these demands is, of course, the de- 
scription of the system by means of appropriate (quan- 
tum) transport equations like the Kadanoff-Baym equa- 
tions. Graphically this is illustrated in fig. ^. In addition 
to the sunset diagram we have also included the Hartree 
diagram there which in a perturbative scheme is the one 
which arises first. The difference to the resummation of 



(4.9) is the fact that the propagators entering into the 
self-energy operators are now also the fully dressed ones. 

Unfortunately, the full quantum transport equations 
are generally hard to solve and thus are not so much 
of practical use. Yet one need not be that pessimistic. 
For weak coupling, i.e. if the damping width is suffi- 
ciently small compared to the quasi-particle energy one 
can take the Markov approximation to obtain standard 
kinetic equations. For the situation of fig. || one gets 



{E p dt - pd s - d s m(x, t)dp) /(£, t; p) = 



(4.10) 



- [fS < (f,i;p, J B p ) (f(x,p,t) + 1) - fE > (x,t;p, J B p )/(x,p,t)] 

Here / denotes the semi-classical non-equilibrium phase- 
space distribution of quasi-particles. m(x, i) denotes 
the sum of the bare and the dynamical (space time de- 
pendent) mass generated by the Hartree term. Within 
this Markovian approximation the fully dressed and re- 
summed propagators are given by [^) 



D Ict / av {x,t-p) 



p2 _ TO 2^ £j _ ReJ]^ £•£,) -j- ip T(x, t;p) ' 

D<(x,t;p) « (-2i) f(x,t;p) 

PoF(x,t;p) 

{p 2 - m 2 (x,t) -ReE(f,t;p)) 2 + p^T 2 (x , t; p) 



(4.11) 
(4.12) 



We emphasize that in (4.12) the instantaneous non- 
equilibrium phase space distribution function f(t) enters 
and not the initial one, h. The dynamically generated 
mass as well as the collisional self energy contribution S 
can thus be evaluated with these propagators. Higher or- 
der terms leading to the pinch singularities are explicitly 
resummed and lead to finite and transparent results. 
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One can now understand the observations made by 
Altherr [^9|. He has found, starting from some non- 
equilibrium distribution h, that higher order diagrams 
contribute to the same order in the coupling constant as 
the lowest order one. Indeed, in his investigation, the 
particular higher order diagrams where not hing but the 
perturbative contributions of the series in (4.4) for the 
dressed or resummed one-particle propagator D < . The 
only difference is that he has employed a 'free' propagator 
modified by some finite width in order that each of the 
terms in the series (4.4) becomes well defined. The reason 
for the higher order diagrams to contribute to the same 
order is that the initial out-of-equilibrium distribution h 
cannot stay constant during the evolution of the system 
as it has to evolve towards the Bose distribution. If n— n b 
is of order o(l), it is obvious that there must exist contri- 
butions which perturbatively attribute to the temporal 
change of the distribution function and contr ibute to the 
same order o(l). In fact, in our prescription ( 4.12 ), h has 
simply be substituted by the actual phase space distri- 
bution /. Then calculating e.g. the tadpole diagram, as 
discussed in the particular case of |19|, one has to stay 
within lowest order in the skeleton expansion, but with 
the fully dressed propagator. 



V. SUMMARY AND CONCLUSIONS 

Our study provides new intuitive insight into non- 
equilibrium quantum field theory and the process of ther- 
malization. In our discussions we have elucidated on 
the stochastic aspects inherent to the (non-) equilibrium 
quantum transport equations, the so called Kadanoff- 
Baym equations. For this we have started with a simple 
model, where we have coupled a free scalar boson quan- 
tum field to an external heat bath with some given tem- 
perature T. We have isolated a term denoted by I which 
solely characterizes the (thermal and quantum) fluctu- 
ations inherent to the underlying transport process. By 
introducing a stochastic generating functional the emerg- 
ing stochastic equations of motion can then be seen as 
generalized (quantum) Langevin processes. Long wave- 
length modes with momenta and energies fc|, lo <C T 
then behave as classical propagating modes for a weakly 
interacting theory. The important observation is that 
there the average occupation number of the soft modes 
becomes large and approaches the classical equipartition 
limit. We hope that our detailed analysis of the real 
time description of the soft modes within the Schwinger- 
Keldysh formalism can attribute to this subject. The 
understanding of the behavior of the soft modes is cru- 
cial e.g. for the issue of baryon number violation in the 
hot electroweak theory || . 

We also have demonstrated how so called pinch sin- 
gularities are regulated within the non-perturbative con- 
text of the thermalization process. These singularities 
do (and have to) appear in the perturbative evaluation 



of higher order diagrams within the CTP description of 
non-equilibrium quantum field theory. Their occurrence 
signals the occurrence of (onshell) damping or dissipa- 
tion. This necessitates in the description of the evolution 
of the system by means of non-perturbative transport 
equations. In the weak coupling regime this corresponds 
to standard kinetic theory. In this case we have given a 
prescription of how the dressed propagators can be ap- 
proximated in a very transparent form. 
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